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The purpose of this note is first of all to explain how intersection numbers of divisors on 
the moduli spaces M. g , n of stable pointed curves can be computed. The Witten conjecture, 
proven by Kontsevich, gives a recipe to compute the intersection numbers of the n basic line 
bundles on A4 g , n • As we will see, the knowledge of these numbers allows one to compute 
all other intersection numbers of divisors as well. That this is possible was pointed out to 
me by Rahul Pandharipande. Earlier, Eduard Looijenga had made a remark that went a 
long way in the same direction. 

After describing the various divisors on A4 3jn , we proceed to discuss the algorithm 
for the computation of their intersection numbers. We have implemented the algorithm 
and we will discuss this implementation as well as the results obtained this way. E.g., we 
computed all intersection numbers on Ai g for g < 6. (Copies of the program and some 
data computed with it are available from the author.) 

A refined version of the algorithm requires to take certain higher codimensional classes, 
introduced by Mumford and Arbarello-Cornalba, into account; it computes all intersection 
numbers of these classes and divisors. Recently, we realized that the Chern classes of the 
Hodge bundle can be taken along as well; hence all intersection numbers of Mumford’s 
tautological classes and divisors can be computed. This has several applications. In §4 
we discuss one application in detail: the calculation of the class of the locus of Jacobians 
in the moduli space of principally polarized abelian varieties of dimension g (projected in 
the tautological ring). This class was classically known for g — 4 and was computed by us 
with ad hoc methods for g — 5; the new method allows in principle to compute the class 
for all g , in practice currently for g < 7. Other applications may be found in the recent 
papers of Graber and Pandharipande [GP] and Kontsevich and Manin [KM], 

§1. Line bundles and divisors on M. g , n . 

For non-negative integers g and n with 2g — 2 + n > 0, denote by A i g , n the moduli 
space of stable n-pointed curves of genus g , over an algebraically closed field k. This 
is the Deligne-Mumford compactification of the moduli space M- g , n of smooth n-pointed 
curves (C; x±,..., x n ) of genus g (with Xi ^ Xj if i ^ j). We will consider certain classes 
in the rational Picard group of -M g ,n ■ First, for 1 < i < n, let ijji denote the first 
Chern class of the line bundle whose fiber at a stable n-pointed curve (C;xi,... ,x n ) is 
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the cotangent space to C at xy; i.e., V’i = °i ( <J i( u; 7 r n+1 )), where 7r n +i : Ad 3 , n _ 1-1 —> Ad 3 , n 
is the morphism obtained by forgetting the (n + l)-st marked point (the universal curve, 
cf. [Kn 1]), u>Tr n+1 is the relative dualizing sheaf, and a±,... ,a n are the natural sections of 
7 r n _|_i (the image of a stable n-pointed curve under a t is the stable (n + l)-pointed curve 
obtained by attaching a 3-pointed rational curve at the i-\ h point and considering the 
remaining 2 points on that curve as the i-th and (n + l)-st point). Next, following [AC 2], 
§1, we define K\ = 7r n+ i,*(/i 2 ), with K = ci(o; 7rn+1 (^” = i Di)), where Di is the divisor 
that is the image of the section ay, for 1 < i < n. Note that it is a consequence of results 
of Harer (cf. [Ha 1], [Ha 2], [AC 1]) that over C the restrictions to Ad s , n of the classes Ki 
and Vh, • • •, V’n generate the rational Picard group of M. g , n ■ 

To get generators for the rational Picard group of Ad g , n • we have to add the fun¬ 
damental classes of the boundary divisors. Exactly when g > 0, there is a boundary 
component whose generic point corresponds to an irreducible singular curve. It is the 
image of Ad 3 _i )n+ 2 under the degree 2 map that identifies the (n + l)-st and (n + 2 )-nd 
point on each curve. Following [AC 2 ], we denote this locus by A ilT and its class in the 
Picard group by hm . (For g — 0 this class is 0 by definition.) 

The other boundary components parametrize reducible singular curves. The generic 
point of such a component corresponds to a curve with two irreducible components C\ 
and C 2 , with genera g\ and (72 satisfying gi -\- g 2 = <7, and labelled by subsets N\ and 
IV 2 of n = { 1 , 2 ,..., rz} satisfying Ni ]j N 2 = n that correspond to the marked points on 
the two components. All partitions of g and n that lead to a stable curve occur; this just 
translates as the condition \N X > 2 when g t = 0. Such a boundary component is the image 
of Ad Sl .| jvi|+i x Ad 92 ,| 7 v 2 |+i under the natural map that identifies the two ‘extra’ points 
and labels the \N t \ remaining points on Ci with the labels from N t . We have chosen to 
denote this boundary component in the case n > 0 by A gij Ni where N t is the subset of n 
containing 1. In the case n = 0 the Ay are empty and we may drop them in the notation; 
note that A S1 = A g2 and that this component is usually denoted as A min (g l 52 ) . 

Although this will play no role in the sequel, we point out that the classes in the 
rational Picard group that we have introduced so far, are independent whenever g > 3. 
For (7 = 2 , there is one relation, originating from the fact that K\ on Ad 2 comes from 
the boundary. For g = 1, both k 1 and the i/y’s come from the boundary; the boundary 
components are independent. For g — 0, the boundary components generate, but are not 
independent; the relations arise from the various projections to Ado ,4 and the equivalence 
of its 3 boundary cycles (cf. [Ke]). Let me add that although the results mentioned in 
this paragraph are no doubt correct, it appears that for some of the statements, references 
containing proofs are not (yet) available. 

There is one other divisor class which will be most useful: Ai, the first Chern class 
of the Hodge bundle. The Hodge bundle is the locally free rank g sheaf (on the moduli 
functor) whose fiber at a curve C is H°(C, wc). So it is 0 in genus 0, while it is a pull-back 
from Ad 1,1 resp. M. g in case g — 1 resp. g > 2 . 

§2. The idea of the algorithm. 

Suppose given a monomial of degree 3g — 3 + n in the divisor classes n\. V’i, • • •, V’n, 
dm and the S t)i .yv ( on Ad g , n ; we want to compute the corresponding intersection number. 
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We will interpret the divisor classes as classes on the moduli functor. In the case of a 
boundary divisor, this means that we divide the usual fundamental class by the order of 
the automorphism group of the generic curve parametrized by the divisor. These divisor 
classes will be denoted by S ..., to distinguish them from the actual boundary divisors A.... 

The case in which the monomial involves only the ipi s is of course covered by the 
Witten conjecture [Wi], proven by Kontsevich [Ko]. As for instance explained in [AC 2], 
this also allows to compute the intersection numbers involving both k\ and the 

It remains to compute the intersection numbers involving a boundary class. Such a 
number may be thought of as the intersection of the remaining classes on the corresponding 
boundary component. A problem with this approach appears to be that most boundary 
components have singularities that are not quotient singularities, which means that one 
cannot properly do intersection theory on them. This problem is easily solved: we have 
seen that each boundary component is the image under a finite map of a moduli space of 
stable pointed curves or a product of two such spaces. (The map almost always has degree 
1; the only exceptions are the degree 2 map from M. g -i , n +2 to A irr and, in the case g even, 
the degree 2 map from M. g / 2,1 x M. g / 2,1 to A g / 2 •) So we wish to pull back the remaining 
divisor classes by means of this map. If we can express the pull-backs in terms of the basic 
classes on the new moduli space(s), we will be done, by induction on the dimension of the 
moduli space. 

So the whole point is to understand the pull-backs of the basic divisor classes from 
Ai gt n to the moduli spaces occurring in its boundary components. 

It is clear that the ipi's pull back to ipiS on the new moduli space(s): to the first n 
on M g -i^n+ 2 , and to the \Ni\ ipi s on M^^n^+i and the |jV 2 | Vh’s 011 -^g 2 ,\N 2 \+i that 
correspond to the points that are not identified in the map to A flij jv 4 . 

As explained in [AC 2], the class pulls back to K\ on Ai g - i,n +2 resp. to the sum of 
the pull-backs of the Kq’s from the two factors on the product x -M P2; |jv 2 |+i . 

Pulling back a boundary divisor different from the one under consideration to the 
new moduli space(s) is not difficult. The main point is that two distinct boundary divisors 
intersect transversally in the universal deformation space (see [DM]). It remains to identify 
the boundary divisors on the new moduli space(s) that arise as the inverse image of the 
intersection of the boundary divisor under consideration with a different one. 

For example, in the case A irr , the pull-back of the class Sh,m to Ai g ~i in +2 is the sum 

4-l,MU{n+l,n+2} + t>h,M > 

with some exceptions: when n — 0 and 2 h = g, the two classes in the sum are equal and 
the pull-back consists of that class just once; when n = 0 otherwise, 6h,® has been denoted 
&g-i-h,{i ,2 } above; when h — 0 or h — g, the first resp. second summand is not defined 
and should be omitted. 

We now consider the pull-backs of boundary divisors to a product 1 ^ 1+1 x 
Aig 2 .\N 2 \+i ■ The pull-back of hm is the sum of the di rr ’s on the two factors. It remains to 
find the pull-backs of the boundary divisors parametrizing reducible curves of a different 
type than the one under consideration. 

We start with the case n — 0. Hence N\ — N -2 — 0. We may assume that g\ < g^ 
and we want to pull back the class 5h , with h < g — h and h ^ (j\ . The general curve 
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in a component of the intersection of the two boundary divisors is a chain consisting of 
3 irreducible components. For the genera of the 3 components, there are a priori the 
following 4 possibilities: 

(1) [h, g i — h, < 72 ] , occurring when gi > h ; 

(2) [h, c /2 — h, gi] , occurring when <72 > h ; 

(3) [gi, h — gi, g — h\ , occurring when h > g 1 ; 

(4) [< 72 , h — r/ 2 , g — h] , occurring when h > g 2 ■ 

Here the second entry refers to the genus of the middle component; note that [a, 6 , c] and 
[c, 6 , a] describe the same type of curves. 

We observe that h < g/2 < g 2 ■ In fact h < g 2 , since equality implies h = g\ , which 
we have excluded. So (4) never occurs, while ( 2 ) always occurs. We conclude that the 
pull-back of Sh equals 

f prfS gi _ hi{ i } + pr* 2 8 g2 _ h , { 1 } in case g x > h ; 

\ P r 2 S g 2 -h,{i} + P r 2 S h- gi ,{i} in case h > g x , 

with one exception: when h = g — h (implying h > g 1 ), the two summands in the second 
line are equal, and the actual pull-back consists of that class just once. 

Now for the case n > 0. We may assume that 1 £ Ni and we want to pull back 
the class S^m ■ different from S gij N x . So 1 G M and (h, M ) ^ (gq, Ni), but we no longer 
have gi < (J 2 or h < g — h. Again there are a priori 4 possibilities for the genera of the 3 
components, where as before the second entry refers to the genus of the middle component: 

(1) [h, gi — h, < 72 ] , occurring when gi > h and M C W ; 

(2) [h, (72 — h, gi] , occurring when g 2 > h and M C N 2 ; 

(3) [gi , h — g \, g — h] , occurring when h > gi and Ni C M ; 

(4) [ 772 , h — 772 , g — h] , occurring when h > g 2 and N 2 C M . 

Observe that (2) never occurs, since 1 e M and 1 ^ Ah. Note also that the other 
possibilities indeed yield stable curves in all cases: in (1) and (3) the necessary condition 
M ^ N\ when h — g 1 is fulfilled, and in (4) the equality M — N 2 never occurs. Finally, 
note that the types (1), (3) and (4) never coincide. 

This means that the pull-back of Sh,M consists of the sum of 0, 1 or 2 of the classes 
from the following list, depending on which conditions are satisfied: 

(1) pr\8i h M (pulled back from A4 fll; jViU{*}) when <71 > h and M C Ah ; 

(3) pr^Sh-g^M-N! u{i} (pulled back from M 92j at 2 u{i}) when h > gx and W C M ; 

(4) pr*5 h _ g2 , M _ at 2 u{*} (pulled back from when h > g 2 and N 2 C M. 

Here we have identified the factors of the product of moduli spaces by means of sets of 
marked points instead of just their number of elements. In light of our convention to label 
a divisor parametrizing reducible curves by the genus of the component containing 1 and 
by the set of marked points on that component, it is natural to give the ‘extra’ point in 
case (3) the label 1, rather than *. 

Finally, we have to deal with self-intersections of boundary divisors: we need to pull 
back the class of a boundary divisor to the corresponding (product of) moduli space(s). It 
is not difficult to deal with this directly, but it is easier to use the fundamental identity 

/■ci = 12Ai — 5 T if) 
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on Jvlg n , where 5 is the sum of the functorial classes of the boundary divisors and is 
the sum of the n i/q’s (see [Co]). Namely, note that every divisor class on M. g , n that we 
have discussed so far, occurs in this identity. So a given boundary divisor class can be 
expressed as a linear combination of other divisor classes, and if we know how to pull back 
the other classes, we will also know how to deal with self-intersections. We have discussed 
the pull-backs of K\ , the ipi s and the other boundary divisor classes above, so we only 
need to determine the pull-back of Ai. Note that the pull-back of the Hodge bundle to 
M g - i,n +2 is an extension of a trivial line bundle by the Hodge bundle in genus g — 1, 
whereas the pull-back of the Hodge bundle to .M^ijvq+i x -M P2 ,|jv 2 |+i is the direct sum 
of the Hodge bundle in genus g\ and the Hodge bundle in genus g 2 ■ (Cf. [Kn 2] or [Co].) 
Hence we find that the pull-back of Ai to A4 g -i t n+2 equals Ai , whereas the pull-back of 
Ai to M gu Ijvp+i X Mg 2y \N 2 \ + i equals pr\\\ +pr^Ai. 

This determines then the pull-back of a boundary divisor class to its corresponding 
(product of) moduli space(s). We find that in the case n — 0 the pull-back of ©r to 
M. g - 1,2 equals 

9~ 2 

—'i/q — '02 + ^irr + ^ , 

h= 1 

whereas in the case n > 0 the pull-back of Ai rr to M. g - i , n +2 equals 

-fpn+i — i>n+2 + ^i rr + (4,MU{n+l> + 4,MU{n+2}) • 

0<h<g-l, lGMCn 
M when h=g — 1 

In the case of a boundary divisor parametrizing reducible curves, an actual self-intersection 
is much rarer. We find that the pull-back of S 91j n 1 to M. gi) N 1 u{*} x -^s 2 ,A r 2 U{i} equals 

- pr \^i + pr 2 S g2 - gi ,{i} when n = 0 and g x < g 2 ; 

-pr*i^{*} ~ pr 2 ipi + pr{5 gi _g 2 , nu{*} when n > 0, N 1 =n and g x > g 2 > 0 ; 

~P r i'4’{*} — P r 2^1 otherwise . 

This finishes the theoretical description of the algorithm. 

§3. Implementation and results. 

We have implemented the algorithm outlined in §2 in Maple 1 . This first of all requires 
an implementation of the algorithm for the computation of the intersection numbers of the 
V’i’s on M. g , n given in [Wi]. For this, we have gratefully used the results of Chris Zaal 
who, using such an implementation, computed a table ([Za]) containing all the intersection 
numbers 

(r dl Td 2 ■■■TdJ = ^1^2 2 

on M g ,n with g < 9 for which all di > 2 (hence n < 24 since ^™ =1 (dj — 1) = 3g — 3). The 
intersection numbers for which a d t equals 0 or 1 can be computed from these by means 
of the string and dilaton equations. 

1 Maple© is a trademark of the University of Waterloo and Waterloo Maple Software. 
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From the intersection numbers of the 'i/’j’s, one can determine the intersection numbers 
of the classes on M. g introduced by Mumford in [Mu], as was briefly explained in [Wi]. 
Arbarello and Cornalba introduced in [AC 2 ] classes Hi on M. g , n that generalize Mumford’s 
classes, and they show that the intersection numbers of the ipi's determine the intersection 
numbers of these h^s as well as the ‘mixed’ intersection numbers of h^s and 'ipi’s. So in 
particular the intersection numbers of the divisor classes Hi and the ipiS are determined. 
It is easy to implement the calculation of these numbers from the intersection numbers of 
the 'i/Cs (especially so with a formula we learned from Dijkgraaf [Dij]), and we have for 
instance calculated the numbers h^ 9 on A4 g for g < 9. 

Naturally the various divisors on M. g ^ n have to be ordered in some consistent way. 
On M. g we start with Hi , followed by hi rr and then the ‘reducible’ boundary divisors, 
ordered by the minimum of the genera of the 2 components, for a total of [g/2] +2 classes. 
(The class Ai was introduced only to deal with self-intersections of boundary divisors and 
is not actually used in the program.) When n > 0, there are (g + l)2 n_1 + 1 classes: 
first V’l , • • •, V’n , «i an d ^irr, then the reducible classes, ordered first by the genus of the 
component that contains the point 1 , then by the number of points on that component, 
and finally by the lexicographic ordering of subsets of n (of equal size and containing 1 ). 
(Recall that in genus 0 the class hi rr is 0; it is included for convenience.) 

In the case of a pull-back to a product of moduli spaces, it is necessary to renumber 
the indices from N\ U {*} as well as those from {1} U IV2 . For this we just use the natural 
ordering of the elements of n , taking * as the (n + l)-st point. 

The implementation of the actual algorithm is now rather straightforward. Given 
a monomial in the divisor classes that contains at least one boundary divisor, we order 
the classes as above, and pull back to the (product of) moduli space(s) corresponding 
to the last occurring divisor. Note that whenever the monomial contains several distinct 
boundary divisors, we have a choice here; while the ordering chosen on M. g is probably 
optimal, this is most likely not so when n > 0. (Only recently have we experimented with 
another ordering, which indeed appeared to be better; the idea is that in the reducible case 
one should try to pull back to two moduli spaces parametrizing curves of approximately 
equal genus and number of points.) 

In case the last divisor is Ai rr , we find a homogeneous polynomial of degree 3g — 4 + n 
in the divisor classes on M . g -\, n +2 • After expanding, it is a sum of monomials (with 
coefficients) in the divisor classes; these are evaluated by means of the (heavily recursive) 
algorithm. 

In the reducible case, we also find a homogeneous polynomial of degree 3g — 4 + n, 
but this time in two sets of variables, the divisors on M. gii \N 1 \+i and those on M- g2 .\N 2 \+i ■ 
Many of the monomials will not have the correct degree 3g\ — 2 + | Ni | in the first set of 
variables and are 0 for trivial reasons. The others automatically have degree 3g^ — 2 + 
in the second set of variables; writing such a monomial as c • Mi ■ M 2 , where c is the 
coefficient of the monomial and Mi is the monic monomial in the i-th set of variables, it 
contributes c • o(Mi) • a(M 2 ), where a(Mj) is the result of applying the algorithm to M,; 
on . 

Using this implementation of the algorithm, we have computed e.g. the 28 intersection 
numbers of K \, d UT and on A4 3 , confirming the results of [Fa 1]. However, because 
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of its heavily recursive character, the algorithm becomes impracticable already in the 
computation of certain intersection numbers on Ad 4 . Most intersection numbers are still 
easy to compute, but especially the numbers n\~ l 8 \ rr with i large take a long time. It 
is quite clear why: firstly, a pull-back to Ad g _ i , n +2 is the ‘worst case’, since the changes 
in genus and in dimension of the moduli space are minimal, while the number of points 
increases by 2, so that the number of divisors increases by a factor of almost 4; secondly, 
as we saw in § 2 , the pull-back of <5i rr to Ad g _i )n+ 2 involves by far the highest number of 
terms. 

Observe however that the class hi rr is a pull-back from Ad 3 resp. Ad q 1 . This first of 
all implies 

5 m+l = 0 , 

where m = max(( 7 , 3g — 3), for all g > 0. A systematic use of this identity already saves 
considerable time. Moreover, any product involving only K\ , the ipi’s and h; rr can be 
pushed down to M. g resp. Ad 1,1 , with the use of the projection formula and the formulas 
in [AC 2 ], This leads to intersection numbers on those spaces of monomials in A irr and 
the higher K t mentioned before. As observed by Arbarello and Cornalba, tlie(ir) K t behave 
very well under pull-back, and it is clear that all intersection numbers involving divisors 
as well as the K t can be computed by means of an algorithm almost identical to the 
one in §2. Even the implementation is easy to adapt. The point is that this greatly 
simplifies the calculation of the numbers involving only K\ , the 'i/y’s and hi rr : in the naive 
implementation, the complexity of a calculation increases at least exponentially with the 
number of points, and one would have to calculate certain numbers on Ado, 2 g in order to 
get all numbers 011 A4 g ; but now many of the hardest numbers can be computed using at 
most 2 -pointed curves at all stages of the computation. 

With these simple changes implemented, the calculation of many more numbers be¬ 
comes practical. We have calculated all intersection numbers of divisors on A4 g for g < 6 
as well as on Adsq and M. 4 , 1 . To obtain the 2 numbers K \and 5^. on Ade , we used the 
relations A} 3 5i rr Ki = 0 and Aj 3 A; 2 rr = 0, consequences of the geometrically obvious relation 
Ai 3 h ir r = 0. 

This section of the paper would hardly be complete without some actual intersection 
numbers. Here are a few: 


O 11 Ad 4 : 
O 11 Ad 5 : 
On Ad 6 : 


9 _ -251987683 x9 _ 1 

irr _ 4320 ’ 1 “ 113400 ’ 

2 _ -1766321028967 12 _ 31 

r_ 6048 ’ 1 _ 680400 ’ 

-32467988437272065977 U5 _ 431 

7257600 ’ 1 “ 481140 ' 


We computed the number Af on Ad 4 in [Fa 2] by a completely different (ad hoc) method. 
Calculating it with the algorithm amounts to calculating all 220 intersection numbers of 
divisors on Ad 4 , so this provides a nice check of the implementation. 
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§4. The class of the locus of Jacobians and other applications. 

The calculation of Af on At 4 in [Fa 2] was used there to obtain the well-known result 
that the class of the locus J\ of Jacobians of curves of genus 4 in the moduli space A 4 of 
principally polarized abelian varieties of dimension 4 equals 8 A 1 . Using the computation 
of A ;} 2 on At 5 as well as some computations in the tautological ring of At 5 as in [Fa 3], 
we could determine the class of in A $, as we will explain in a moment. Recently we 
realized that Mumford’s formula [Mu] for the Chern character of the Hodge bundle on 
M g , together with an algorithm similar to the one discussed in § 2 , enable one to compute 
the class of J g in A g for all g. at least in principle. In practice, we have carried this out 
for g < 7. These results will be discussed here as well. 

First we recall the set-up and explain what we mean by “the class of the locus of 
Jacobians in the moduli space of principally polarized abelian varieties”: this will be the 
class in the tautological ring of A g , the Q-subalgebra of the cohomology ring of a toroidal 
compactification A g of A g generated by the Chern classes A t of the Hodge bundle E on 
A g . From [Mu], §5 we know that the relation 

(1 — Ai + A 2 — A 3 + ... + (—l) 3 A g ) (1 + Ai + A 2 + A 3 + ... + A g ) — 1 

holds; equivalently, c/i 2 fc(E) = 0 for all k > 1. The tautological ring is in fact the quotient 
of Q[Ai,..., A g ] by the ideal generated by the homogeneous components of the relation 
above and is thus a complete intersection ring. A detailed description of the tautological 
ring may be found in [vdG]. In particular, the relation with the cohomology ring of the 
compact dual of the Siegel upper half space via Hirzebruch’s proportionality principle is 
explained there. This includes the fundamental identity 

9 9 l r l 

iR.=n^ 

i —1 2=1 

that enables one to compute intersection numbers in the tautological ring of A g . 

Denoting by t : M g —>■ A g the extended Torelli morphism and its image by J g , we are 
after the functorial class [ J 9 \q of the locus of (generalized) Jacobians, which is one half its 
usual fundamental class. In other words, we wish to determine 1, since a generic curve 
of genus at least 3 has 110 non-trivial automorphisms, while the generic p.p.a.v. and the 
generic Jacobian of dimension at least 3 have two automorphisms. 

It is important to point out that this is not what we actually compute. We do not 
know whether the class of the locus of Jacobians actually lies in the tautological ring (our 
feeling is that this is not the case for g large enough, but we cannot even think of a method 
to decide this). Instead, we compute the projection of this class in the tautological ring; 
this is well-defined by the perfect pairing in the tautological ring and the cohomology ring 
of A g . (In other words, we compute the class modulo a class X that pairs zero with all 
classes in the tautological ring of the complementary dimension.) 

The method of computation of this class is the following. It is a class in the tautological 
ring of A g of dimension 3g — 3, hence of codimension c = ( 9 ~ 2 ). Write it as a linear 
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combination with unknown coefficients a % of the elements s, of a basis of the degree-c part 
of the tautological ring: 

[J g \ Q — 1^*1 = a lSl + (12S2 + ... + cikSk (+-X")- 

(A natural basis is for instance the collection of square-free monomials of degree c in 
Ai,..., X g .) To compute the coefficients Oj, we have to evaluate the k monomials A, t of 
a basis of the degree-( 3 g — 3 ) part of the tautological ring on this class. The evaluation 
of expressions A t s : j in the tautological ring of A g uses the relations between the Aj and 
the proportionality relation stated above. So the expressions Aj(aiSi + ... + cikSk ) yield k 
linearly independent rational linear combinations of the unknowns a,. The values of these 
expressions can be determined as -jt*(f*Aj) = (Af* 1) ■ Aj, if we know how to evaluate t* Aj 
on M g . 

The simplest non-trivial example is g = 4. Here c = 1 , a basis in codimension 1 is Ai, 
a basis in codimension 9 is Af (or any 11011-zero monomial in the Aj), so to compute the 
class of J,\ in A4 we only need A ! ( on M4. We have seen above that this can be evaluated 
e.g. by means of the implementation of the algorithm; we find the well-known result that 

[JAq — 8A1. 

The situation for g = 5 is more interesting. Here c = 3 and a basis in codimension 3 is 
given by A1A2 (= ^Af) and A3 . We need to evaluate 2 independent monomials of degree 12 
in the Aj on M 5 . The algorithm will naturally yield only the number A{ 2 (whose value we 
gave at the end of § 3 ). However, we can use the simple observation that the class \ g \ g -i 
vanishes on the boundary M g — M g (see [Fa 3 ] or [Fa 4 ]). As a corollary, the numbers 
A1A2A4A5 and A3A4A5 on M 5 satisfy the same relation as the classes A1A2 and A3 in the 
1 -dimensional degree -3 part of the tautological ring R*(M5) of M 5. This relation was 
worked out in [Fa 3 ]: IOA3 = 3A1A2. A quick calculation in the tautological ring of A$ 
shows that this implies that the class of the locus of Jacobians satisfies 

[Jb ]q — a (3A1A2 — 2A3). 

Another such calculation, using X\ 2 = 6g j^ 00 on M 5 , shows then that a = 24 , hence 

[J§\q = 72A1A2 — 48A3. 

For higher genus, we find with this method only some of the coefficients, not all of 
them. We start with a general formula: 

Conjecture 1. In the basis of monic square-free monomials in A1 ,X fJ of degree 2 ), 
the coefficient Ci^,...,g-2, of A1A2 • • • A 3 _3 in the (projected) class [J 9 \q equals 

1 3-2 2 

1 TT_ t _. 

2(7-2 f \ (2i + l)|S 2 j| 

So it equals 1 , 8, 72, 384, 768 resp. for g = 3 ,..., 7, while it is not an integer for any larger 
value of g. 

The conjecture is true for g < 15 , as it is a consequence of a conjectural formula in [Fa 3 ] 
for the number A^_ x on M g that is proven for g < 15 . To derive it from that formula, 
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simply note that = 2 A 3 _ 2 A g _iA g and that A 1 A 2 • • • A fl _3 is the unique monomial in 

the basis that pairs non-zero with A s _ 2 A s _iA g . 

Note that we could have used this formula instead of that for X \ 2 to compute the 
class of [Jfj]q ■ 

For g — 6 we Hnd the coefficient of A2A4 using the relation between A1A3 and A4 in 

R*(M 6 ): 

C 2A = - 3 ( 71 , 2,3 = -1152. 


The knowledge of Ap (see §3) gives a non-trivial relation between the remaining coefficients 
of A 1 A 5 and Aq: 


Cq + 16 ( 71,5 — 


7336704 

691 


A new ingredient will be required to solve for these coefficients. It is provided by 
Mumford’s formula [Mu] for the Chern character of the Hodge bundle on A4 g : 


ch{ E) = g + ^ 

i=1 


R‘ 2 i 

m 


1 9 ~ l 

«2t-l + { K V~ 2 

Z h =0 


— K 2 i ~ 3 K 2 + ... + Kf- 2 ) 


(Note that we use a different convention for the Bernoulli numbers: B 2 — B 4 = -^ etc.) 
Here Iq : A4 g -y j2 —»• Aj rr C A4 g and i\ x : M.h,i x ■M. g -h,i —^ A/j C M. g are the natural 
maps, and K, is the first Chern class of the relative cotangent line bundle at the z-tli point. 

The formula is ideally suited for a recursive computation of the intersection numbers 
of the K t and the ch 2 j-i(E). Namely, suppose given a monomial in those classes of degree 
3g — 3 + n on M. g , n ■ If only K t occur, we can proceed as explained in §2, using the Witten 
conjecture (Kontsevich’s theorem). In any case, as the Hodge bundle is a pull-back from 
A4 g resp. Mi,i , we can push down the expression to M. g resp. M. 1,1 and obtain a sum of 
similar expressions. In genus 1, we only need the well-known equalities 

chi(E) = Ai = «i = Vh = (ti) = ^ • 


So assume the genus is at least 2 and the monomial on A4 g contains at least one ch 2 j-i(E). 
Take the highest odd Chern character component that occurs, and expand it using Mum- 
ford’s formula. In the first term, that ch 2 k~i is replaced by a n 2 k~i (up to a factor), so it 
is determined inductively. The other terms involve expressions in the Kj pushed forward 
via the maps iy, . The point is that these can be written as push-forwards from M . g -\, 2 
resp. M.h .,1 x A / tg-h,i of intersection numbers of the classes K, , Kj and Chern character 
components of the Hodge bundles in genus g — 1 resp. genera h and g — h. This is clear 
from the Arbarello-Cornalba formulas for and the fact that Zq( 1E) is the extension of 
a trivial line bundle by the Hodge bundle in genus g — 1, while for h positive i£(E) is the 
direct sum of the Hodge bundles in genera h and g — h. 

After expanding and omitting the terms that are 0 for dimension reasons, we find an 
expression in intersection numbers of the classes just mentioned on spaces of 1 - or 2 -pointed 
curves. These can be pushed down again to intersection numbers of K t and ch 2 j-i(K) on 
M.h (with h < g) or A4 1,1 . By induction on the genus, these numbers are known. 
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Having discussed the implementation of the divisor algorithm in some detail, we con¬ 
tent ourselves with saying that the implementation of the new algorithm proceeds along 
similar lines and is considerably easier. 

Because the A i can be expressed in the chj (E), this means that all the intersection 
numbers of the A i on M. g can be computed recursively. By the discussion above, it follows 
that for all g the projection in the tautological ring of A g of the class of the locus of 
Jacobians can be computed, at least in principle. 

Currently, we have carried this out for g < 7. For g = 6 . only one more relation was 
required. Either one of the two relations following from 


A2A3A4A6 — 


1697 


2988969984000 ’ 


A1A2A3A4A5 


150719 

15692092416000 


suffices to solve for Cq and (the first, relation involves Cps only). The result is 

248064 


474048 

[Jg]q = 384 A 1 A 2 A 3 — II 52 A 2 A 4 -|-—A 1 A 5 


— 2 3 ( A 1 A 2 A 3 — 3 A 2 A 4 + 


691 
2469 

-Ai Ar 

1382 


691 


- Af 


646 

691 


Af 


It may be worthwhile to point out the relation 15Cq + 28(7^5 = 2 9 3 3 . 
In genus 7 the result is: 


\Ji]q — 768 A 1 A 2 A 3 A 4 — 69 I 2 A 2 A 3 A 5 -|-—-A 1 A 4 A 5 -|-—-A 1 A 3 A 6 


691 


691 


8842752 968832 3276672 

-A4A6 H-——A3A7-——A1A2A7 . 


691 


691 


691 


As stated, this result is not very pretty; perhaps the class looks better in a different basis. 
We would like to point out though that the class can be computed with any choice of 
7 independent monomials of degree 18 in the A i ; once this is accomplished, the values 
of all other such monomials can be determined by means of an easy calculation in the 
tautological ring of A- . In particular, one can choose ‘easy’ monomials to compute the 
class, and get the ‘hard’ ones for free. In this way we computed for instance A j 8 on A 4 -: 


A 


18 

1 


32017001 
638512875 ' 


Finally, we have another general formula: 

Conjecture 2. In the basis of monic square-free monomials in A1 ..... X g of degree , 
the coefficient C2,3,..., g -4, g -2 of A2A3 • • ■\ g -4\ g -2 in the (projected) class [, J 9 \q equals 


( - 2 ) _ 29 - z \ r = — FT 2 _ 1 _ TT 4 

V 12 J 1 ’ 2 ”"’ 9-3 12 f \ (2i + l)\B 2 i\ Ag-A f = \ (2i + l)\B 2i 
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So it equals -48, -1152, -6912 resp. for g — 5, 6 , 7, while it is not an integer for any larger 
value of g. 

Again, the conjecture is true for g < 15; it is a consequence of Conjecture 1 and of the 
conjectural formula n\X g - 3 — g(2g — 2)A 3 _ 2 in R 9 ~ 2 {M . g ) that can be proven for g < 15 
using the results of [Fa 3] (note that the only 2 monomials in the basis that pair non-zero 
with Ai X g — 3 X g — i X g are A 2 A 3 • • • A g _ 4 A g _ 2 and AiA 2 .. . A g _ 3 ). 
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Finally, some remarks about possible uses of the implemented algorithms. The al¬ 
gorithm for computing intersection numbers involving Chern classes of the Hodge bundle 
(described in §4) has already found applications in the recent work of Graber and Pandhari- 
pande [GP] and Kontsevich and Manin [KM]. For moduli spaces M. g , n of small dimension 
and with not too many divisors, one can use the divisor algorithm to determine the part 
of the cohomology ring generated by divisors. This makes e.g. an intersection calculation 
on All .4 as in [Ge] relatively painless. As in that paper, there may be applications to com¬ 
puting Gromov-Witten invariants. In fact, the algorithms involving the K t and the c/i 2 j_i 
allow to do calculations that include those classes; this covers e.g. the Chow ring of AI 3 . It 
is even possible to include arbitrary boundary strata as module generators, by pulling back 
all other classes to the corresponding product of moduli spaces via a sequence of maps, 
each identifying a single pair of points. Writing algorithms that can handle intersections 
of arbitrary boundary strata will be considerably more difficult, however. 

The various algorithms, and some tables of intersection numbers computed with it, 
are available from the author by e-mail, although only the initial version of the divisor 
algorithm is currently available in a user-friendly format. 
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